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ABSTRACT 

The statistical expectation values of the temperature fluctuations of cosmic 
microwave background (CMB) are assumed to be preserved under rotations of the 
sky. This assumption of statistical isotropy (SI) of the CMB anisotropy should 
be observationally verified since detection of violation of SI could have profound 
implications for cosmology. We propose a set of measures, k £ (£ = 1, 2, 3, . . .) for 
detecting violation of statistical isotropy in an observed CMB anisotropy sky map 
indicated by non zero K£. We define an estimator for the Kg spectrum and ana- 
lytically compute its cosmic bias and cosmic variance. The results match those 
obtained by measuring K£ using simulated sky maps. Non-zero (bias corrected) 
Kg larger than the SI cosmic variance will imply violation of SI. The SI mea- 
sure proposed in this paper is an appropriate statistics to investigate preliminary 
indication of SI violation in the recently released WMAP data. 

Subject headings: cosmic microwave background - cosmology: observations 

The Cosmic Microwave Background (CMB) anisotropy is a very powerful observational 
probe of cosmology. In standard cosmology, CMB anisotropy is expected to be statisti- 
cally isotropic, i.e., statistical expectation values of the temperature fluctuations AT(g) 
are preserved under rotations of the sky. In particular, the angular correlation function 
C(q, q') = (AT(q)AT(q')) is rotationally invariant for Gaussian fields. In spherical harmonic 
space, where AT(q) = Y.i m a imYimi.q) this translates to a diagonal (ai m a*, m ,) = Ci5i V 5 mm > 
where Ci, the widely used angular power spectrum of CMB anisotropy, is a complete de- 
scription of (Gaussian) CMB anisotropy. Hence, it is important to be able to determine if 
the observed CMB sky is a realization of a statistically isotropic process, or not 1 . 



1 Statistical isotropy of CMB anisotropy and its measurement has been discussed in literature ear- 
lier (Ferreira & Magueijo 1997; Bunn & Scott 2000). 
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We propose a set of measures, K£ {t = 1,2,3,.. .) which for non-zero values indicate and 
quantify statistical isotropy violation in a CMB map. A null detection of Kg will be a di- 
rect confirmation of the assumed statistical isotropy of the CMB sky. It will also justify model 
comparison based on the angular power spectrum C\ alone (Bond,Pogosyan & Souradeep 1998,2000). 
The detection of statistical isotropy (SI) violation can have exciting and far-reaching impli- 
cation for cosmology. In particular, SI violation in the CMB anisotropy is the most generic 
signature of non-trivial geometrical and topological structure of space on ultra-large-scales. 
Non-trivial cosmic topology is a theoretically well motivated possibility that is being observa- 
tionally probed on the largest scales only recently (Ellis 1971; Lachieze-Rey & Luminet 1995; 
Starkman 1998; Levin 2002). 

For statistically isotropic CMB sky, the correlation function 



C(n u n 2 ) = C(ni ■ h 2 ) = — J dTZCiTZh^ Tln 2 ), (1) 

where TZn denotes the direction obtained under the action of a rotation 1Z on n, and dIZ is 
a volume element of the three-dimensional rotation group. The invariance of the underlying 
statistics under rotation allows the estimation of C(n\ ■ 77,2) using the average of the temper- 
ature product AT{n)AT{n') between all pairs of pixels with the angular separation 9. In 
particular, for CMB temperature map AT(^) defined on a discrete set of points on celestial 
sphere (pixels) (i — 1, . . . , N p ) 



C{9) = Xr(q t )&r(q J )5(cos9 - q, ■ q 3 ) , (2) 

is an estimator of the correlation function C{9) of an underlying SI statistics 2 . 

In the absence of statistical isotropy, C(q, q') is estimated by a single product AT(q)AT(q') 
and hence is poorly determined from a single realization. Although it is not possible to es- 
timate each element of the full correlation function C(q,q'), some measures of statistical 
anisotropy of the CMB map can be estimated through suitably weighted angular averages of 
AT(q)AT(q'). The angular averaging procedure should be such that the measure involves 
averaging over sufficient number of independent 'measurements', but should ensure that the 
averaging does not erase all the signature of statistical anisotropy (as would happen in eq. (1) 
or eq. (2)). Another important desirable property is that the measures be independent of the 



2 This simplified description does not include optimal weights to account for observational issues, such as 
instrument noise and non uniform coverage. However, this is well studied in literature and we evade these 
to keep our presentation clear. 
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overall orientation of the sky. Based on these considerations, we propose a set of measures 
Kg of statistical isotropy violation given by 

1 2 



K t = dtt dtt' 



(21 + 1) 
8tt 2 



-J dn X e(n)C(nq,nq') 



(3) 



where C(lZq, IZq') is the two point correlation between IZq and IZq' obtained by rotating 
q and q' by an element 7Z of the rotation group. The measures Kg, involve angular average 
of the correlation weighed by the characteristic function of the rotation group XeiJZ) = 
D e MM (1Z) where D e MM , are the Wigner D- functions (Varshalovich et al. 1988). When 
7Z is expressed as rotation by an angle uo ( where < uo < it) about an axis r(0,<3>), the 
characteristic function Xiffi) = Xe( u ) = s ^ n [(2^ + l)u;/2]/sin [uj/2] is completely determined 
by uo and the volume element of the three-dimensional rotation group is given by dlZ = 
4 sin 2 uo/2duo smOdOd®. Using the identity / dTZ'xe(TZ')xe(TZR') = xe(R>), expression (3) 
can be simplified to 

= JdnJdn>c(q, q ')J dn X i(n)C(n q ,n q :) (4) 

containing only one integral over the rotation group. For a statistically isotropic model 
C(1Zqi, TZq 2 ) = C(cji, q 2 ) is invariant under rotation, and eq. (4) gives Kg = K°5g , due to 
the orthonormality of Xei^)- Hence, k,£ defined in eq. (3) is a measure of statistical isotropy. 

The measure Kg has clear interpretation in harmonic space. The two point correlation 
C(q, q') can be expanded in terms of the orthonormal set of bipolar spherical harmonics as 

C(q, q') = A uf{ Y M) ® Y v {q')} m , (5) 

U'iM 

where Ajjf are the coefficients of the expansion. These coefficients are related to 'angular 
momentum' sum over the covariances (ai m a*, m ,) as 

A uf = ^( a /m. a /*'m/) ( — 1)™ £lml'-m> ; (6) 
mm' 

where £fmi'm' are Clebsch-Gordan coefficients. The bipolar functions transform just like 
ordinary spherical harmonic function YgM under rotation (Varshalovich et al. 1988). Substi- 
tuting the expansion eq. (5) into eq. (3) we can show that 

«< = £W>0, (7) 

U'M 

is positive semidefinite and can be expressed in the form 

2£ 



K t 



8tt 2 



-fdfcXt(K) Yl ( a l™ a Vm>)( a lma*> m ') n i ( 8 ) 
J Iml'm' 
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where (. . ) n is computed in a frame rotated by 1Z. When SI holds (a« m ^* m ') = Ci^iv^mm 1 ', 
implying = (— 1)'C;(2Z + 1) 1/2 5 H / <5mo- The coefficients Aff represent the statistically 
isotropic part of a general correlation function. They also represent the statistically isotropic 
part of any arbitrary correlation function. The coefficients Af^f are inverse-transform of the 
two point correlation 

A™ = JdnJ dQ'C(h, ti) {Y h (n) ® Y l2 (n')}* m . (9) 

The symmetry C(h,h') = C{h',n) implies 

A m = ( _ 1)( * 1+ , 2 -^m A m = A m ^ fc = 0> ^ 2> . . . . (10 ) 

Recently, the Wilkinson Microwave Anisotropy Probe (WMAP) has provided high res- 
olution (almost) full sky maps of CMB anisotropy (Bennet et al. 2003) from which k,£ can 
be measured. Given a single independent CMB map, AT(g) we need to look for violation 
of statistical isotropy. Formally, the estimation procedure involves averaging the product 
of temperature at pairs of pixels obtained by rotating a given pair of pixels by an angle ui 
around a sufficiently large sample of rotation axes. The integral in the braces in eq. (3) is 
estimated by summing up the terms for different values of u weighed by the characteristic 
function. We can define an estimator for Kg. as 



(21+1) 
8tt 2 



J2 AT(g,) ATiqj) ^ X e(w m ) ^ AT(ft mn %)AT(ft mn ^) , (11) 



m=l 



n=l 



where as described below 23 ^ = (kf) accounts for the 'cosmic bias' for the biased estimator 
kf . As with the sky, the rotation group is also discretized as 7Z mn where m — 1, . . . ,N W 
is an index of equally spaced intervals in rotation angle w and n — 1, . . . , N r indexes a set 
of equally spaced directions in the sky. While we have also implemented this real space 
computation, practically, we find it faster to estimate Ke in the harmonic space by taking 
advantage of fast methods of spherical harmonic transform of the map. In harmonic space, 
we first define an unbiased estimator for the bipolar harmonic coefficients based on eq. (6) 
and then estimate Kg using eq. (7) 
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(12) 
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Assuming Gaussian statistics of the temperature fluctuations, the cosmic bias is given 
by (Hajian & Souradeep 2003) 

(kf ) = Yl [« mi a «im' 1 )«m 2 %^) + (at mi a l2m > 2 )(al m2 a hm{ )] 

h,h mi,m'j rri2,m 2 

x X^<r m <r m , , (id 

x ^M 1 mi(2m 2 N 1 m' 1 !2m^ ■ 

M 

Given a single CMB sky- map, the individual elements of the (ai m aj, m ,) covariance are 
poorly determined. So we can correct for the bias that arises from the SI part of corre- 
lation function where 

e+h 

fB t = (kf) si = (2£+ 1) E C *A i 1 + i-V'Shh] ■ (14) 

h h=\e-h\ 

Hence, for SI correlation, the estimator ke is unbiased, i.e., (ke) = 0. 

Assuming Gaussian CMB anisotropy, the cosmic variance of the estimators A\jf and kg 
can be obtained analytically for full sky maps. The cosmic variance of the bipolar coefficients 



im 

) (ai 2m2 a l2m ' 2 ) + (a hmi ai 2m > 2 ) (ai 

m\,m' 1 m2,m' 2 

which, for SI correlation, further simplifies to 



vlMtl) = c h c l2 [i + (-iY5 hh ] Y: (-ir^c&^C^- (16) 

mi ,m2 

Note that for li = l 2 the cosmic variance is zero for odd £ due to eq. (10) arising from 
symmetry of C(q, q'). 

A similar but more tedious computation of 105 terms of the 8-point correlation function 
yields an analytic expression for the cosmic variance of kg (Hajian & Souradeep 2003). For 
SI correlation, the cosmic variance for t > is given by 



l:2l>£ 



2 ^rrf + ( ~ 1} ' (2£ + 1} + (1 + 2 (- 1 ^ F » 
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n 2 

e+h 

E * 

_«2 = K-il| 

+ i6(-D' E i f± i f E c ?a. < 17 > 

h h I 

„.l,„ ro _ \^ \^ \^ r lM r lM r tM' r lM' /io\ 

Wlieie -T;^ — ^ /-^ ^h-rn 1 h-rn 3 ^hm 2 hrn i ^hm, i hm 1 { ^h-rn2,h-rn 2 - 

m\mi=—l\ m3tr»4=- I3 M,M'=—£ 

Numerically, it is advantageous to rewrite F^, in a series involving 9-j symbols. The expres- 
sions for variance and bias are valid for full sky CMB maps. For observed maps one has to 
contend to incomplete or non uniform sky coverage. In such cases one would estimate the 
cosmic bias and variance from averaging over many independent realizations simulated CMB 
sky from the same underlying correlation function. Fig. 1 shows the measurement of K£ in a 
SI model with flat band power spectrum. The bias and variance is estimated from making 
measurements on 50 independent random full-sky maps using the HEALPix 3 . The cosmic 
bias and variance obtained from these realizations match the analytical results. Just as in 
the case of cosmic bias, the cosmic variance of K£ at odd multipoles is smaller. The figure 
clearly shows that the envelope of cosmic variance for odd and even mulitpole converge at 
large £. For a constant /(/ + 1)C; angular power spectrum the <j si (kh) falls off approximately 
as \ji at large £. (The absence of dipole and monopole in the maps affects K£ for £ < 4 
leading to the apparent rise in cosmic variance at £ < 4 seen in Fig 1.) 

The bias and cosmic variance depends on the total SI angular power spectrum of the 
signal and noise Ci = Cf + . Hence, where possible, prior knowledge of the expected 
signal should be used to construct multipole space windows to weigh down the contribution 
from the region of multipole space where SI violation is not expected, e.g., the generic 
breakdown of statistical isotropy due to cosmic topology. The underlying correlation patterns 
in the CMB anisotropy in a multiply connected universe is related to the the symmetry of 
Dirichlet domain (Wolf 1994; Vinberg 1993). In a companion paper, we study the K£ signal 
expected in flat, toroidal models of the universe and connect the spectrum to the principle 
directions in the Dirichlet domain (Hajian & Souradeep 2003). SI violation arising from 
cosmic topology is usually limited to low multipoles. A wise detection strategy would be to 
smooth CMB maps to the low angular resolution. When searching for specific form of SI 
violation, linear combinations of K£ can be used to optimize signal to noise. Before ascribing 
the detected breakdown of statistical anisotropy to cosmological or astrophysical effects, one 



3 Publicly available at http://www.eso.org/science/healpix/ 



h h=\£-h\ h 1 
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must carefully account for and model into the SI simulations other mundane sources of SI 
violation in real data, such as, incomplete and non-uniform sky coverage, beam anisotropy, 
foreground residuals and statistically anisotropic noise. 

In summary, the Kn statistics quantifies breakdown of SI into a set of numbers that can 
be measured from the single CMB sky available. The Ke spectrum can be measured very fast 
even for high resolution CMB maps. The statistics has very clear interpretation as quadratic 
combinations of off-diagonal correlations between a\ m coefficients. Signal SI violation is re- 
lated to underlying correlation patterns. The angular scale on which the off-diagonal correla- 
tions (patterns) occur is reflected in the Ke spectrum. As a tool for detecting cosmic topology 
(more generally, cosmic structures on ultra-large scales), the K£ spectrum has the advantage 
of being independent of overall orientation of the correlation pattern. This is particularly 
suited for search for cosmic topology since the signal is independent of the orientation of the 
DD. (However, orientation information is available in the A\^f .) The recent all sky CMB 
map from WMAP is an ideal data set where one can measure statistical isotropy. Interest- 
ingly, there are hints of SI violation in the low multipole of WMAP (Tegmark et al. 2003; 
de Oliveira-Costa et al. 2003; Eriksen et al. 2003). Hence is of great interest to make a care- 
ful estimation of SI violation in the WMAP data via Ki spectrum. This work is in progress 
and results will be reported elsewhere (Hajian et al. 2003). This approach complement di- 
rect search for signature of cosmic topology (Cornish, Spergel & Starkman 1998). 

TS acknowledges enlightening discussions with Larry Weaver, Kansas State University, 
at the start of this work. TS also benefited from discussions with J. R. Bond and D. Pogosyan 
on cosmic topology and related issues. 
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Fig. 1. — The figure shows the bias corrected 'measurement' of K£ of a SI CMB sky with a flat 
band power spectrum smoothed by a Gaussian beam (l(l+l)Ci = exp(— Z 2 /18 2 ) . The cosmic 
error, ct(k^), obtained using 50 independent realizations of CMB (full) sky map match the 
analytic results shown by lower dotted curve with stars. The upper dotted curves separately 
outline the cosmic error envelope for odd multipoles (filled triangles) and for even multipoles 
(empty triangles) to explicitly highlight their convergence. Violation of SI will be indicated 
by non-zero Ke measured in an observed CMB map in excess of a(ne) given by the C\ of the 
map. The lower dashed curve (filled squares) shows the cosmic error for ideal unit flat band 
power spectrum (1(1 + 1)C/ = 1) with no beam smoothing. The curve falls off roughly at l/£ 
at large £. 



